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In contrast to recent studies [Rabl, Phys. Rev. Lett. 107,063601 (2011); Nunnenkamp etal, Phys. Rev. Lett. 
107, 063602 (2011)] on photon blockade that prevents subsequent photons from resonantly entering the cavity 
in optomechanical systems, we study the photon induced tunneling that increases the probability of admitting 
subsequent photons in those systems. In particular, we analytically and numerically show how two- or three- 
photon tunneling can occur by avoiding single-photon blockade. Our study provides another way on photon 
control using a single mechanical resonator in optomechanical systems. 

PACS numbers: 42.50.Wk, 07.10.Cm, 42.50.Ar, 42.65.-k 



I. INTRODUCTION 

Optomechanical systems (HQ] have attracted extensive at- 
tention in the past few years because of its potential appli- 
cation in high-precision measurements and quantum informa- 
tion processing. To realize these benefits, the mechanical res- 
onator has to be in its ground state; and the optomechanical 
radiation-pressure interaction strength should be bigger than 
the decay rates of the cavity field and the mechanical oscil- 
lator. The ground state cooling of the mechanical resonator 
has been experimentally studied in optomechanical systems 
(e.g., in Refs. Ikl UToll ). Although the latter condition is not 
easy to be satisfied in standard optomechanical systems, how- 
ever the experimentalists have obtained an effectively strong 
coupling, by applying a classical driving field to the cavity 
mode of the optomechanical systems, which has led to obser- 
vations of normal-mode splitting (e.g., in Refs. lflTl[l2"tO and 
optomechanically-induced transparency (e.g., in Refs. lfl3l - 
However, the resulting effective coupling resembles two 
linearly coupled harmonic oscillators lUPl in their optical and 
mechanical position quadratures, and the coupling strength is 
proportional to the square root of the mean cavity photon num- 
ber, thus it does not really describe the nonlinear effect at the 
single-photon level. 

Recent two theoretical proposals lfl6i IT7I] showed that the 
single-photon effect or photon blockade IU8I1 can occur when 
the optomechanical systems are approaching the strong cou- 
pling at the single-photon level H 1 9142211 . This is because the 
frequency of the cavity field is parametrically modulated by 
the motion of the mechanical resonator in the optomechanical 
systems, this modulation results in the photon-photon inter- 
action lf23Tr28ll . When the strong optomechanical radiation- 
pressure interaction makes the nonlinear photon-photon cou- 
pling strength bigger than the decay rate of the cavity field, 
then the photons can prevent the subsequent photons from 
resonantly entering the cavity. The photon blockade phe- 
nomenon has been demonstrated experimentally in cavity 



QED systems for microwave Ir29[|30| and op tical ll3lll32ll pho- 
tons. Meanwhile, the experiments I3U 13211 also observed the 
photon induced tunneling, that is, the probability of admit- 
ting subsequent photons is increased when there is one pho- 
ton inside the cavity. Moreover, absorption and emission of 
resonant photons in pairs have been observed in a single atom 
strongly coupled to an optical cav ity IT331 l34ll . 

Motivated by these works IU6U341I an d also in contrast to 
the studies for photon blockade IU61I17I1 . here we investigate 
photon-induced tunneling in the optomechanical systems. In 
Sec. II, we will describe the theoretical model of the system, 
and introduce the master equation, which is used to study sta- 
tistical properties of the cavity field. The effect of the mechan- 
ical resonator on the mean photon numbers is discussed via 
the master equation when the optomechanical system reaches 
the steady state. In Sec. Ill, we will study the two-photon tun- 
neling via the normalized second-order coherence function of 
the cavity field by numerically solving master equation. In 
Sec. IV, we further use the three-photon tunneling as an ex- 
ample to show how the multi-photon tunneling can happen. 
Finally, a summary is given in Sec. V. 



n. THEORETICAL MODEL 

We study an optomechanical system (see, review in 
Ref. (U]) in which a cavity field is coupled to a mechanical 
resonator through the radiation pressure. The Hamiltonian can 
be written as 



H = tuoQa) a + huj m b^b + hGa) a Uy + 6) 



(1) 
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where a (b) and cv (&') are the annihilation and creation 
operators of the cavity field (mechanical resonator with the 
fundamental mode) with the frequency ojq (w m ). The pa- 
rameter G describes the coupling between the cavity field 
and the mechanical oscillator. To obtain the eigenvalues of 
the Hamiltonian in Eq. (Q~|), let us apply a unitary transform 
U = exp [—Ga^a (p* — b) /oJ m ] to Eq. fl}, and then we ob- 
tain an effective Hamiltonian H' = U HU^ with 

/ Q2 \ Q2 

H' = huj m b^b + h I ujq ] a) a — h a)a)aa, (2) 
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which has eigenstates \n, m) and corresponding eigenvalues 

E n .m = h (nuiQ - n 2 —) + hmuj m , (3) 

V w m y 

where \n,m) = U\n,m), \n,m) represents that there are n 
photons and m phonons in the optomechanical system, re- 
spectively. Because a unitary transform does not change the 
eigenvalues of the Hamiltonian, thus the eigenvalues of the 
Hamiltonian in Eq. (Q~|) are the same as in Eq. (01. 

We now assume that the cavity field is coupled to a weak 
probe field with the frequency u c and the coupling strength 
e c . In the rotating reference frame with the unitary operator 
R(t) — exp[iui c (r at], the Hamiltonian in Eq. (0 becomes 

H = hAa^a + hui m tfb + %Ga)a (tf + b) + ihe c (a f - a) ,(4) 

where A = uq — uj c is the detuning between the cavity field 
and the probe field. When the environmental effect is taken 
into account, the dynamical evolution of the reduced density 
operator p for the cavity field and the mechanical resonator 
can be described via the master equation ll35ll 



dp 

dt 



1 

ih 



H,p 



- (2apa) 



a) ap — pa) a) 



(2bptf - tfbp - ptfb) 



+7 m n m (bpb f + tfpb - tfbp - pbtf) 



(5) 



with the decay rates 7 and j m of the cavity field and me- 
chanical mode. n m is the mean phonon number of the me- 
chanical resonator given by the Bose-Einstein statistics n m = 
[exp(huj m / ksT) — with the Boltzmann constant ks and 
the environmental temperature T of the mechanical resonator. 
Here, the frequency of the cavity field is assumed to be high 
enough such that the thermal photon effect can be neglected. 
In the basis of the photon and phonon number states \n, m), 
the formal solution of the reduced density operator p(t), gov- 
erned by the master equation in Eq. (|5j, can be given by 



P(t) 



n,m,n f ,m' 



n.m:n ,m 



(£)|n,ra)(n', m'\ 



(6) 



If all elements p n ,m;n' ,m'(t) in Eq. (0 are given, then any 
physical quantity of the optomechanical system can be calcu- 
lated. For example, the mean photon number inside the cavity 
can be calculated as 
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FIG. 1 : (Color online) The mean photon number (n) inside the cav- 
ity as a function of the normalized detuning A / Ao in (a) for differ- 
ent decay rates of the mechanical mode at the given environmental 
temperature T = 1 pK and (b) for different temperatures and given 
decay rate j m /2Tr = 0.01 MHz. Other parameters in (a) and (b) are 
£ c /2tt = 0.01 MHz, G/2-k = 2.5 MHz, 7/2tt = 0.1 MHz, and 
cu m /2TT = 10 MHz. 



is maximum at (A/Arj) = 1 which corresponds to the fre- 
quency uj c = ujo — (G 2 /w m ) of the probe field. As shown in 
Eq. (O, this is because that the transition frequency of the cav- 
ity field from the one-photon state 1} to the vacuum (ground) 
state |0) is shifted to luq — (G 2 /u! m ) when the mechanical res- 
onator is coupled to the cavity field. Additionally, Fig. [Tta) 
and (b) show that when the decay rate j m or the environmen- 
tal temperature T of the mechanical resonator increases, the 
full width at half maximum in the curve of the mean photon 
number (n) becomes broad for other given parameters. This 
means that the decay and environmental temeperature of the 
mechanical resonator significantly affect the lift time of the 
cavity photons. 



III. TWO-PHOTON TUNNELING 

The tunneling, blockade and other statistical properties of 
the photons can be characterized by the normalized nth-order 
correlation function at the zero time delay 



0<">(O) = 



(a tn a") = TrQt"q") 
(at a )" = Tr(pata)" ' 



(8) 



In our study here, the reduced density operator p of the cavity 
field and the mechanical resonator can be obtained by solving 
the master equation in Eq. (0. 

Let us now study the simplest photon tunneling in pairs. 
Similar to the photon blockade, this phenomenon can be char- 
acterized by the normalized second-order correlation function 



(n) = (0+0) = Tv\p{t)a^a] = V ¥ „ m . n m (t) . (7) 9 W (0) at the zero time delay. When gW (0) is less than one 



Using Eqs. (O and (0 and in the steady-state condition, the 
mean photon number (n) is plotted as a function of the nor- 
malized detuning A/Ao with Ao = G 2 /u m for: (i) differ- 
ent decay rates of the mechanical mode, e.g., j m /2n = 0.01 
MHz, 7 m /27r = 1 MHz, and 7 m /2?r = 2 MHz in Fig. [J 
(a) and (ii) different environmental temperatures of the me- 
chanical mode, e.g., T = 1 fjK, T = 5 mK, and T = 10 
mK in Fig.[T](b). We find that the mean photon number (n) 



the photons inside the cavity block subsequent photons from 
resonantly entering the cavity, and single-photons come out 
of the cavity. However, when (0) is bigger than one, pho- 
tons inside the cavity enhance the resonantly entering prob- 
ability of subsequent photons, this photon induced tunneling 
can make photons come out of the cavity in pairs under the 
certain condition. In Refs. ll33ll36ll . the second-order differen- 
tial correlation function at the zero time delay 



(0) = (atV) - {a) a y = [g^ (0) - 1] 



(9) 
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FIG. 2: (Color online) g (2) (0) and C (2) (0) as functions of the nor- 
malized detuning A / Ao and the coupling strength G for T = 1 /iK 
in (a) and (b), respectively, g' 2 ' (0) as a function of A/Ao for dif- 
ferent temperatures T at G/2-k = 2.5 MHz in (c). C (2) (0) as a 
function of A/Ao for different temperatures T at G/2n — 0.75 
MHz in (d). Other parameters are e c /27r = 0.01 MHz, 7/27r = 0.1 
MHz, uj m /2ix = 10 MHz, and 7 m /2?r = 0.01 MHz. 



is introduced to characterize the probability of creating pho- 
ton pairs simultaneously in the cavity, where (n) is the mean 
photon number inside the cavity. 

To explore two-photon tunneling and compare it with the 
photon blockade, we study g^ 2 '(0) and C^ 2 '(0) in the steady- 
state. In Figs. Sa) and (b), 5 (2) (0) and C<- 2 \0) versus the de- 
tuning A — ujq — uj c and the coupling strength G are plotted at 
the environmental temperatures T = 1 /iK of the mechanical 
resonator. In each figure, there are two curves which corre- 
spond to single-photon (green dashed curve for A = Ao) and 
two-photon (black dashed curve A = 2A ) resonant excita- 
tions from the ground state to the first excited and the second 
excited states of the cavity field, respectively. Figs. |2ja) and 
(b) clearly shows g {2) {0) < 1 and C^{0) < for the single- 
photon resonant excitation from the ground state to the first 
excited state when the coupling strength G is bigger than the 
decay rate 7 of the cavity field, that is, the single-photon phe- 
nomena or photon blockade occurs as shown in Refs. |[T6[[T7ll . 
However, if the frequency of the probe field equals to the half 
of the transition frequency from the ground state to the second 
excited state, i.e., A = 2Ao, then the single-photon transition 
from the ground state to the first excited state is suppressed, 
but the second photon can enter the cavity and makes resonant 
transition from the ground state to the second excited state to- 
gether with the first photon. This photon-induced tunneling 
means that the cavity can absorb photons in pairs simultane- 
ously, therefore we can find <?' 2 '(0) > 1 and C' 2 '(0) > in 
Figs. Ha) and (b). 

We further showed the effect of environmental temperature 
on the photon tunneling via ,g' 2 '(0) and C^ 2 '(0) in Figs. [2jc) 
and (d). We find that the low temperature can make photon- 
induced tunneling and photon blockade easy to be observed, 
because with the decrease of the temperature, <p 2 '(0) and 
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FIG. 3: (Color online) (a) Energy level diagram for the optomechan- 
ical system. Here, the horizontal long black lines linked by the ver- 
tical black lines denote the energy levels of the cavity field with fre- 
quency cjo- The parameter ui c is the frequency of the driven field. 
A = wo — w c . The horizontal bright purple lines denote the energy 
levels (|n, m)) of the optomechanical system with Ao = G 2 /w m . 
The horizontal long blue lines linked by the vertical blue line with 
two arrows denotes the transition of |1, m) —¥ |2, in), (b) <7' 2 '(0) is 
plotted as functions of G/u m for A = A , T = 1 /iK. (c) ff (2) (0) 
is plotted as functions of G/u) m for A = 0, T — 1 fiK (black) and 2 
mK (red). The other parameters are the same as in Fig. |2jc). 



C^(0) decrease under the single-photon (A = Ao) reso- 
nant driving, and increase under the two-photon (A = 2Ao) 
resonant driving. As shown in Fig. [2 a), there is an additional 
resonant peak for g^ 2 '(0) at A = 3Ao, which corresponds to 
the transition between two eigenstates |1, fh) and |2, m). We 
find that the height of the resonant peak at A = 3 Ao increases 
when the temperature becomes higher, qualitatively, this is be- 
cause the environmental temperature can change the popula- 
tion distribution, especially it can enhance population distri- 
bution in higher energy levels, and then the transition from 
\l,fh) to j 2, fh) is also enhanced. All these result in the in- 
crease of (0) at A = 3 Ao when the temperature becomes 
higher. 

The effect of the phonon states on the photon blockade has 
already been studied in Refs. lfl6ll27ll . It was found that un- 
der the photon blockade condition A = Ao, the two-photon 
resonant transition between |0, 0} and |2, fh) can also occur 
if Ao = mu! m /2. Thus there are resonant peaks located at 
G/uj m = y/m/2 for g^ (0) as shown in Fig.Qb) for differ- 
ent phonon states |m) (e.g., to = 1,2,3,4). The phonon states 
can also affect the resonant condition of the photon-induced 
tunneling. As shown in Fig. Oc), under the resonant condi- 
tion A = 0, there are five resonant peaks for <J^ 2 '(Q) 3> 1 
with G < oj m . As shown Fig. Oa), these five peaks cor- 
respond to two types of the resonant conditions: (i) if Ao 
satisfies the condition 4Ao = mu) m , then two-photon reso- 
nant excitation between |0, 0} and |2, fh) can occur with the 
peaks located at G/oj m = Wm/i, m = 1,2,3; (ii) if Ao 



4 



satisfies the condition 3Ao = muj m , then resonant excita- 
tion from |1, 0} to |2, to) is allowed with the peaks located at 
G/uj m = y/m/Z, m = 1,2. We also find that the transition 
of 1 1 , 0) — > 1 2, to) is enhanced with the increase of the tem- 
perature, so the resonant peak at G/ui m — \J to/3, to = 1, 2 
increases as the temperature becomes higher. 

To further analyze the conditions of the photon blockade 
and the photon induced tunneling, we now study the relation 
between g^ (0) and the probabilities P(n) of n-photon distri- 
bution. Because the probabilities P (n) corresponding to the 
state in Eq. © are P (n) = ^2 m Pn, m -n,m(t), then we have 
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FIG. 4: (Color online) (a) g (2) (0) and 2P(2)/P(1) 2 , (b) p (3) (0), 
6P(3)/P(1) 3 , ^ 3) (0) and 3P(3)/[P(1)P(2)] plotted as functions 
of the normalized detuning A/A for G = 2.5 MHz and T = 1 /iK. 
The other parameters are the same as in Fig. |2jc). 



IV. MULTI-PHOTON TUNNELING 



In the limit of the weak probe field, e.g., e c = O.Ty, P (n) ^> 
P (n + 1) and in the steady-state, g^ 2 '(0) can be given ap- 
proximately as 



,(2) 



(0) 



2P(2) 

WW 



■i(A 



7 + i2(A-2A ) 



(ID 



here we have assumed that the phonon is in its ground state 
to obtain second approximated expression. This assumption 
has also been made in following derivations in Eq. (15[ and 
Eq. (TH . Eq. ( fTTT i means that g^ 2 \0) is proportional to the 
ratio between the probability to prepare a two-photon state 
and that to prepare two single-photon states independently. In 
Fig.Ua), g( 2 > (0) and its approximated expression in Eq. (flU 
versus the normalized detuning A are numerically simulated 
in the steady-state. We find that the approximated expres- 
sion fits well with the exact solution of g( 2 \0) in the limit 
of the weak probe field. Therefore, g& (0) > 1(5^(0) < 1) 
means that the probability to excite two-photon state is bigger 
(smaller) than that to excite two single-photon states indepen- 
dently, and then photon-induced tunneling (photon blockade) 
happens. To further explore the relation between the height 
of Fig.|2jc) (and also Fig. Ufa)) and system parameters, from 
Eq. dTTb . we can approximately obtain 



In principle, the multi-photon tunneling (more than two) 
can be studied via Eq. ([8]). However, for the sake of simplic- 
ity, we just use the normalized third-order correlation function 
t/ 3 )(0) as an example to show multi-photon tunneling. Using 
the same discussions as forEq. ( fTQb and Eq. ( fTTT i. g^ (0) can 
be approximately written as 



9 i3) (0) 



Tr (pa^a 3 ) 6P (3) 



[Tr(pata)]3 [ P 



(15) 



6[i+z(A-A )] ; 



[7 + i2 (A - 2A )] [§ 7 -M3 (A - 3A )] 



in the steady-state case with the weak probe field. It is clear 
that g^ (0) is proportional to the ratio between the probability 
to create a three-photon state and that to create three single- 
photons independently. If we introduce a quantity 



9? (0) - 9 



(3) 



(0) 



Tr (pa 



gW (0) Tr (pat a) Tr (pa^a 2 ) 



(16) 



to characterize the ratio between the normalized third-order 
and second-order correlation functions, then in the steady- 
state case, we can approximately obtain 



3 (2) (0) 



7 2 +4Ag 

7 2 ' 



(12) 



for A = 2Ao, which corresponds to the height of the peak 
in Fig. |4](a) and Fig. Etc) at A = 2A . Thus the height of 
the peak in Fig. Ufa) is approximately determined by the ratio 
between the decay rate 7 of the cavity field and the parameter 
Aq. We can approximately obtain another formula 



9 {2) (0) 



7 



V 



4A§' 



(13) 



for A = An, which corresponds to the deep in Fig. |4](a) and 
Fig.|2jc) at A = An. It is clear that the conditions to observe 
photon blockade or photon induced tunneling respectively are 



4A 2 > 7 2 



or A > 7/2. 
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(17) 



which is proportional to the ratio between the probability to 
create a three-photon state and the joint probability to create a 
single-photon and a two-photon state independently, which is 
also a possible way to create three photons. 

In Fig. H|b), 5 ^(0) and its approximated expression 
6P(3)/[P(1)] 3 as well as g^(0) and its approximated 
expression 3P(3)/[P(1)P(2)] are plotted as functions of 
A/ An. We find that the approximation in both Eq. (TT~5T > and 
Eq. (fTTI i is valid for the weak probe field. Under the condition 
we can further find: 



that 4Ag > 7 2 
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(18) 



5 



for the resonant condition A = Ao; and 

5 (3) ( o)«^l, 5 f(0)«l, (19) 

for the resonant condition A = 2A ; as well as 

,<•> 9 f(0),^, ,20, 

for the resonant condition A = 3Ao. Eqs. (TT9ll20l ) show 
that the heights of the resonant peaks are approximately deter- 
mined by the ratio of Ao and 7 as for two-photon tunneling. 

Fig. gib) and Eqs. ( 118H20t show that: (i) for the single- 
photon resonant excitation (A = Ao), g^{0) < 1 and 

92 (0) < 1' tn i s tells us that the photon inside the cavity 
is antibunching or the subsequent photons will be blocked by 
the photons inside the cavity, (ii) For the three-photon reso- 
nant excitation (A = 3A ), g (3) (0) > 1 and g { 2 3) (0) > 1, that 
is, the photon inside the cavity is bunching or the cavity can 
absorb three photons simultaneously, (iii) Remarkably, for 
the two-photon resonant excitation (A = 2Ao), <?^(0) > 1 
which means that the probability of three-photon absorption 
is bigger than that of three single-photon absorption indepen- 
dently. We also find gf^ (0) ss 1 at the point A = 2A . 
This tells us an expected result that the joint probability of 
the photon absorption by the cavity in pairs after or before 
single-photon absorption is approximately equal to that of 
three-photon absorption by the cavity. Therefore, the neces- 
sary condition for absorbing three photons by the cavity si- 
multaneously (or three-photon tunneling) is A = 3Ao- In this 
condition, g (3) (0) > 1 and g£' (0) > 1, the probability of 



generating three photons simultaneously is bigger than those 
of generating three single-photons independently and gener- 
ating a single-photon after or before generating a photon pair. 



V. CONCLUSIONS 

In summary, we have studied the photon induced tunnel- 
ing and compare it with the photon blockade in optomechani- 
cal systems with the single-photon strong coupling condition. 
Our study shows that the cavity field can exhibit photon anti- 
bunching when the probe field is resonant with the transition 
from the ground state to the first excited state of the cavity 
field. However, the two-photon tunneling occurs when the fre- 
quency of the driving field equals to the half of the transition 
frequency from the ground state to the second excited state 
of the cavity field. Moreover, we also find that three-photon 
tunneling occurs when the frequency of the probe field satis- 
fies the condition of the three-photon resonant excitation from 
the ground state to the third excited state. Our studies also 
show that the phonon states greatly affect the multi-photon 
resonance in the strong optomechanical coupling at the single- 
photon level. These studies can be easily generalized to the 
multi-photon tunneling case. Our results show that it is pos- 
sible to observe photon induced tunneling in the optomechan- 
ical systems under the single-photon strong coupling condi- 
tion. 
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